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We study restoration of spontaneously broken symmetry at nonzero temperature in the framework 
of the 0(2) model by using polar coordinates. We apply the CJT formalism to calculate the 
masses and the condensate in the double-bubble approximation, both with and without a term that 
explicitly breaks the 0(2) symmetry. We find that, in the case with explicitly broken symmetry, 
the mass of the angular degree of freedom becomes tachyonic above a temperature of about 300 
^ ' MeV. Taking the term that explicitly breaks the symmetry to be infinitesimally small, we find that 

■ the Goldstone theorem is respected below the critical temperature. However, this limit cannot be 

performed for temperatures above the phase transition. We find that, no matter whether we break 
the symmetry explicitly or not, there is no region of temperature in which the radial and the angular 
degree of freedom become degenerate in mass. These results hold also when the mass of the radial 
' mode is sent to infinity. 

PACS numbers: 11.10.Wx,12.39.Fc,11.30.Rd 

I. INTRODUCTION 

O 
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Ch . The 0{N) model is of great irnportance for the theory of critical phenomena []| and has been intensely studied in 
different theoretical frameworks [2t9|. Using a mass term with a 'wrong' sign, the 0{N) symmetry is spontaneously 
broken to 0{N — 1) in the vacuum. At sufficiently high temperature T, the 0{N) symmetry is restored. The 
QQ ■ case = 4 has been profusely investigated because of its application to the chiral phase transition in quantum 
0^ . chromodynamics (QCD) l3"15|. QCD has a chiral SU{Nf)B, x SU{Nf)L symmetry which is spontaneously broken to 
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SU{Nf)ii+L in the vacuum. The group 0(4) is locally isomorphic to SU{2) x SU{2), and 0(3) is locally isomorphic 
to SU{2). Therefore, using the universality hypothesis, the 0(4) model can be considered as an effective theory for 
the restoration of chiral symmetry in QCD with two quark flavors. 

Most studies of symmetry restoration within the 0{N) model 10-1^. have been performed in Cartesian 



coordinates, — ((/)i, 02, • ■ • , (f'N)- Then, when the symmetry is spontaneously broken, the first coordinate is usually 
selected to assume a non- vanishing vacuum expectation value tp, (f)i —> ip + (jji, and the remaining N — 1 coordinates 
, are taken to be the N — I Goldstone modes arising from spontaneous symmetry breaking. However, in the absence of 
M ' explicit symmetry breaking terms, the effective potential is of the well-known 'Mexican hat'-type, with any state on 
the circle |(/)| = = const, being energetically degenerate with the vacuum state (p. Therefore, the Goldstone modes 
actually correspond to angular variables when moving along the circle \(f)\ = const., and the massive degree of freedom 
corresponds to the radial degree of freedom when trying to 'climb' the rim of the hat. 

This observation constitutes the main motivation for the present paper where we study symmetry restoration within 
the 0{N) model in polar coordinates. For the sake of simplicity, we shall concentrate on the case N ~ 2, because 
the transition from Cartesian to polar coordinates is particularly simple in this case. Nevertheless, the conceptual 
issues are the same in this case as for N = 4. In a certain sense, the case N = 2 corresponds to QCD with one flavor 
in the absence of the axial anomaly. For Nf = 1, the (most simple) anomaly term would introduce a linear term 
in the massive degree of freedom, corresponding to an explicit symmetry breaking term. This is different from the 
two-flavor case, i.e., A'^ = 4, which corresponds to the case of maximal anomaly or maximal U{\)a breaking. Due 
to the similarities with the chiral symmetry of QCD and its breaking and restoration at nonzero T, in the following 
we shall also refer to the transition in the 0{2) model as 'chiral transition'. The order parameter for the transition, 
(/?, will be called 'chiral condensate', and the two degrees of freedom (the 'chiral partners') are referred to as the a 



meson and the pion. The so-called 'chiral limit' is the case where explicit symmetry breaking terms are sent to zero, 
in which case the pion becomes a true Goldstone boson. 

Another motivation for our work is that, if we send the mass of the radial degree of freedom (the a meson) to 



infinity, we just obtain the loading-order Lagrangian of chiral perturbation theory (xPT) |16l - ll8j . with the Goldstone 
boson (the pion) as angular degree of freedom. In our simplified framework of the 0{2) model, it is then possible to 
perform a test of the validity of angular variables at nonzero T. 

The transition from Cartesian to polar coordinates represents a change of the representation of the same theory. 
Obviously, physical quantities should be independent from the adopted representation. This fact is ensured by the 
so-called S-matrix equivalence theorem [19[ which states that the elements of the S'-matrix do not change when 
performing a transformation of the fields. However, for actual calculations one must use a certain approximation 
scheme. Then, the results obtained in one representation are not necessarily equal to those obtained in another 



representation 2C|. In this work we shall explicitly show that, by using the CJT formalism in the double-bubble 
21[, quantities computed in polar coordinates are actually quite different from those evaluated in the 



approximation 
standard Cartesian coordinates. 

Our results are the following: (i) for explicit symmetry breaking, the pion mass becomes tachyonic at high T, 
signalizing a breakdown of the CJT formalism, (ii) When the explicit symmetry breaking term is sent to zero (which, 
in analogy to the case = 4, we refer to as the chiral limit), the Goldstone boson (pion) remains massless at each 
T in polar coordinates, thus satisfying the Goldstone theorem. This is contrary to Cartesian coordinates where a 
nonzero value of the pion mass is obtained as soon as the temperature is switched on, see e.g. Ref. [12[. However, due 
to singular terms ~ in the gap equations, and because the order parameter ip = above the phase transition, the 
model becomes ill-defined, (iii) Both with and without explicit symmetry breaking, there is no region of temperature 
in which the chiral partners become degenerate in mass, (iv) Similar results hold also in the nonlinear limit, i.e., when 
the radial excitation becomes infinitely heavy. 

The reason why polar coordinates are problematic at high T can be traced back to the decreasing value of the 
chiral condensate ip: in fact, in polar coordinates there are interaction terms proportional to inverse powers of ip, 
which render the application of the CJT formalism (or any other resummation scheme) problematic when ip is too 
small. In order to circumvent this problem one can perform a slightly different transformation to polar coordinates, 
in which the polar coordinates are not defined with respect to the origin of the Cartesian coordinates 22|. In this 
way a smooth limit from Cartesian to polar coordinates is realized, in which all the results of Cartesian coordinates 
can be reobtained. Interestingly, it is possible to investigate these issues also in the very simple situation of a free 
Lagrangian, see Sec. II VI for details. In the Cartesian representation the results for thermodynamical quantities, such 
as the pressure, are exact in this case. The deviations of the results in the polar representation from the Cartesian 
one explicitly show the limitations of polar coordinates. 

A further issue of polar coordinates is the fact that the Jacobian associated with the field transformation is not 
unity: an additional term emerges in the transformation of the interaction measure, which is potentially relevant in 
the context of quantum field theory. In order to investigate this point we study a different field representation in 
terms of polar variables, in which the Jacobian is indeed unity. The qualitative picture does not change and the same 
conceptual issue of exploding interaction terms in the high-temperature region exists also in this case. 

The paper is organized as follows: In Sec. |ll]we write the 0{2) model in terms of polar coordinates. In Sec. IIIII we 
apply the CJT formalism and present the numerical results. The simple case of a free Lagrangian is discussed in Sec. 
IIVI In Sec. |V] the alternative representation with unit Jacobian is discussed. Finally, we outline our conclusions and 
give an outlook in Sec. IVII Our units arc fi. = c = fcs = 1; the metric tensor is g^jy = diag (+, —,—,—). 
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II. THE 0(2) MODEL IN POLAR COORDINATES 



A. Tree level, zero temperature 

The 0(2) model in Cartesian coordinates (f> = {(f>i,(j)2), including an explicit symmetry breaking term H, is 
described by the Lagrangian 

1 ^ _> 772^ X -* -> 

Heart = ^d^,^ • 9^0 + — • (/) - -((/) • + H<t>^ . (1) 

As usual, we consider the shift (^i — )■ </)i + where (/) is a constant. At zero temperature the minimization of the 
potential V[(\>\ = (/),(/)2 = 0) leads to the minimum 4> = satisfying the following gap equation: 

mV - 2A(/33 + H = . (2) 

When H = the symmetry of the Lagrangian ([T]) under 0(2) transformations is exact. For > and H ^ 
the global minimum is realized for cj) = ip ^ 0, i.e., the ground state breaks the 0(2) symmetry spontaneously. The 
vacuum expectation value (p is referred to as 'chiral condensate' in our model. In the vacuum, the numerical value for 
ip is chosen to be the pion decay constant, ip = f.^ ~ 92.4 MeV. When H ^ an additional explicit breaking of chiral 
symmetry is realized. 

By shifting the fields around their values at the global minimum, one obtains the zero-temperature tree-level masses 
nil and m2 as 

ml = -m2 + 6A/2 , mj = -m^ + 2\fl . (3) 

It is clear that m2 — > for i/ — > 0, i.e., this particle represents the Goldstone boson emerging from the spontaneous 
breaking of chiral symmetry. 

We now introduce polar coordinates (a, tt) through the transformation 



leading to 



TT 7T 

bi = a cos — , 4>2 = cr sin — , (4) 



Cpoi = Wad^a + —d^nd'^Tr + '-^g' - '-a^ + Ha cos'- . (5) 



Just as above, one shifts the field a as a ^ a + ((>. At zero temperature the minimization of the potential leads to 
the same Eq. ([2]) for the minimum (f> = ip. Also the zero-temperature tree-level masses and to^ coincide with the 
expressions of Eq. ^ : 

2 2 2 , n \ -[2 2 2 ^ /r\ 

= nil = + "A/^ ; = ttIj = — . (o) 
In order to extract m^, we have expanded the cosine in Eq. ([S])- 

B. Mathematical issues using polar coordinates: The Jacobian and the integration intervals 

Denoting Tr/(j) = 9 and taking into account that the Jacobian of the transformation in Eq. (j4]) is a, the partition 
function can be rewritten as 



l/T \ OO 277 

Z = i> 'D(f>iD(j)2 exp I — / dr I dxCE.cart = f a f Vd exp 



dT j dxCE.poi I , (7) 



J V ° 



V 
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where periodic boundary conditions are understood: 0i(t, x) = (f)i{T + 1/T, f), etc. The sufSx E means that the 
Lagrangians are considered in EucUdean space. 

The r.h.s. of Eq. ([7]) describes the partition function in polar coordinates. Due to the fact that the Jacobian is 
not unity and that both fields a and 6 do not vary between (— cx),C)o), the question arises if we can apply the usual 
Feynman rules to the Lagrangian Cpoi- 

One can rewrite the contribution of the Jacobian as 



oo / l/T 

Va (7 — > <j) Va exp 
y 



-A^ / dT j dx\n^\ (8) 



where A is a UV cutoff. It is evident that this term induces a divergent contribution to the effective action requiring 
regularization. In the framework of dimensional regularization the contribution of the Jacobian vanishes in virtue of 



Veltman's rule [23[, see also the explicit perturbative analyses performed in Refs. [2^l2J]. In this work, however, we 
follow a different way. In the present section we neglect the contribution arising from the Jacobian. This is at first 
sight not legitimate because we do not work in dimensional regularization, but with different regularization schemes 
(trivial and counter-term regularizations, see Appendix 1X1 and Sec. IIIIF)) . In Sec. |V]we introduce polar coordinates 
in a slightly modified manner which corresponds to a unitary Jacobian. We shall then a posteriori verify that the 
features described in the present section are qualitatively unchanged also when including the effects of the Jacobian. 
We now turn to the extension of the range of integration over the fields a and tt. The possibility to extend the 

2tz oo 

angular integration interval, J T>6 — > J DO, originates from the periodicity of the integrand. This point is subtle 

-oo 

since it is in general not possible to split a path integral over the interval I = IiU I2 into the sum of two path integrals 
over the intervals Ji and I2, respectively. For potentials of the form U{(j) and for 27r-periodic potentials [/{a, 6) (as 
in the present case) one can show that extending the range of integration simply yields a countably infinite overall 



25| 



constant which can be absorbed into a normalization constant 

The extension in the a direction from (0, 00) to (—00,00) can be achieved with the help of a modified Heaviside 
step function defined in such a way that its contribution vanishes in dimensional regularization [2^ . Note that in this 
work either the trivial or the counter-term regularizations arc used. Although wc do not expect a qualitative change 
of our results, we plan to perform a study using dimensional regularization for the counter terms. 



C. Shift of the potential 

In this section we show how to circumvent the problems arising from the fact that polar coordinates are ill-defined 
at the origin. Inspired by Ref. [22[, we first shift the potential along the 0i-axis by an arbitrary amount v > 0. In 
this way the global minimum realized at the critical temperature (and above) is not located at = 0, but at cj) = v. 
After the shift, the Lagrangian ([T|) reads 



= Mid^q^i + M2d^^2 + ^ (01 - vy + —4>l 



A r 



2 2 

When performing the transformation to polar coordinates, Eq. we obtain 



Cv = Cpoi — rn^v a cos ( 



1 



\v{v - 2cr cos 9) (2(7^ + - 2av cos ( 



Hi 



-v) 



Hv 



(9) 



(10) 



Note that obviously Cv=o — Cpoi from Eq. (O, thus the study of Sec. Ill Al can be regarded as a special case of this 
more general treatment. After performing a Taylor expansion of the trigonometric functions about 9 = and shifting 
(T — !■ (T + 0, we can easily determine the additional tree-level contributions to the masses and the interaction vertices. 
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III. RESULTS AT NONZERO T 



111 this section we present the results at nonzero temperature for the model described by the Lagrangian £„ for 
different cases. 



A. H^O,v^O 



In this case the system is described by the Lagrangian ([5]) in Sec. Ill Al The effective potential in the CJT formalism 
reads 



H 



where U denotes the classical potential and the inverse tree-level propagators read 

At nonzero T the condensate and the masses become T-dcpcndcnt functions 

and the dressed propagators are given by 

G^i(fc) = ~k^ + m2 , G-i(fc) = -Z2fc2 + A/2 , 



(11) 



(12) 



(13) 



(14) 



where Z is a wave-function rcnormalization factor for the pioii. 

The term V2 in Eq. ([TT|) is the contribution of all 2PI vacuum graphs, <j> denotes the connected 1-point function in 
the presence of a source, and G denotes the full connected 2-point function in the presence of the source. In general, 
V2 consists of infinitely many diagrams, which prohibits an explicit calculation of V2. In practice, one therefore has 
to restrict oneself to certain classes of diagrams. We shall use the so-called Hartree approximation where only the 
double-bubble diagrams arc taken into account: 



V2 = -A 



G,(fc) 



H 



G.(fc) 



1 

202 



By extremizing the effective potential in Eq. (jlip we obtain the following equations: 



2,. o^,.3 , JL ^G.(fc) + ^ 



mV - 2A^3 -^H = Q\ip / GJk) - „ , 

' 2(p'' 





2 








j G^{k) 


1 






j PG^{k) 


.k 




-9 




.k 



Ml ^ -m^ + 6A^2 Qx I Gaik) ~ ^ I k^G^ik) 



k 



= 1 



1 
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Gaik) 



(15) 
(16) 
(17) 
(18) 



where (f) = (f denotes, in general, an extremum. At T = the masses coincide with their tree-level values in Eq. ([3]) 
and if — /t^. Furthermore, 



(19) 
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FIG. 1: Upper row {H 7^ 0), from left to right: masses, extrema of the effective potential, pressure (dashed: Cartesian, solid 
lines: polar, 1; = 0). Lower row {H = 0), from left to right: masses, extrema of the effective potential, pressure (dashed: 
Cartesian, solid lines: polar, v = 0). The dotted and dash-dotted lines in the middle panel of the upper row correspond to the 
condensate in the polar case with v = O.S/tt and v — 0.7325/,r, respectively. 



A = 



4(^2 



(20) 



= . (21) 

The following numerical values are used at zero temperature: mo- = 600 MeV, m^r = 139.5 MeV, and = 92.4 MeV. 

In this section wc present the results using the so-called trivial regularization: the vacuum part of the integrals is 
simply set to zero, for details sec Appendix [A] For a discussion of the alternative counter-term regularization scheme 
we refer to Sec. IIIIFI 

The numerical results for the masses, the condensate, and pressure can be found in the upper row of Fig. [1] The 
results of polar coordinates are also compared to the standard results obtained in Cartesian coordinates. One observes 
that in general the results in polar coordinates differ substantially from those in Cartesian coordinates, except for 
the pressure at low temperatures, where the agreement between solid and dashed curves is seen to be quite good. 
Moreover, although the condensate decreases sharply at a temperature ^ 300 MeV, indicating the onset of chiral 
symmetry restoration, the pion mass does not become degenerate with the a mass for high T. On the contrary, above 
a temperature T^ax — 300.5 McV (where the curves in Fig. [1] terminate) it becomes tachyonic, signalizing a breakdown 
of the model for large T. This fact shows a limitation of the application of angular variables at nonzero temperature. 

The effective potential is shown in Fig. [2] for different temperatures. There is a region around the origin where the 
effective potential is not defined since no real-valued solutions exist. Due to the singular terms with inverse powers of 
(j) there is no extremum at the origin. At low temperature there is only a global minimum at a large value = $ > 0. 
At a certain temperature a second minimum and a maximum (both at smaller values of 0) occur. At T* ~ 279.6 MeV 
both arc at the same level, indicating a first-order phase transition. Above this temperature, the global minimum 
<i> moves closer and closer to the origin, but never becomes zero. Above T„iax — 300.5 MeV no real solutions to the 
system of equations exist at the global minimum. 

The reason why no real solution exists above Tmax is due to the fact that the pion mass becomes imaginary, and 
therefore the system is not stable. In order to see this in more detail, we show the function A/7r(0) at T Tmax 
in Fig. |3] (solid line). We observe that M^((/)) has an imaginary part below a value (j>turn — 15 MeV. At (ptum the 
imaginary part vanishes and becomes a positive definite, real- (but multi-) valued function of At T — Tmax, 
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the value at (ptum coincides with at the global minimum of the effective potential. The solid line above (ptum is 
the solution for M.^{(l)) which has to be used to plot the effective potential. Below T^^ax the point 4>turn is located 
such that the extremum Lp > (pturm i-C, real solutions at the global minimum $ exist. With increasing temperature 
$ approaches ipturn and finally hits the ill-defined region for Tmax- 
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FIG. 2: H ^0,v - 
T* = 279.6 MeV. 



0. Effective potential minus its value at its global minimum <j) ~ ^ for temperatures T = 278. 9f MeV and 




10 20 30 40 50 

[MeV] 



-1.826X 10*' 

-1.828x10' 

,"-1.83x10' 
> 

S;1.832x 10' 
'*-1.834x lo' 
-1.836x10' 
-1.838x10' 




15 20 25 30 35 40 

<p [MeV] 



FIG. 3: H ^ 0, V — 0, T — 300.5 MeV ~ Tmax- Upper figure: Pion mass as function of the variable (p. Solid line: MTr{<j)) when 
only the real solution exists, cf> > 4>turn- Dashed line: ReM.^ (</>), dotted line: ImM^ (</;). Lower figure: effective potential as 
function of the variable (j). The latter can be plotted only for (j) > cf>turn, i.e., where real- valued solution for the pion mass exist. 
For 4> < 4>turn also the effective potential becomes a complex function. 
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FIG. 4: H = 0,v ^0. Effective potential minus its value at its global minimum for T = (dotted), T = 0.1 MeV, T = 0.999/^, 
T* = 2.787/^ (solid lines) and T = 2.87/^ (dotdashed) . 



B. H = 0, v = 

In this section we study the case where the explicit symmetry breaking parameter is set to zero, H = 0. Equations 
igi-dTll) simpHfy to 



if - 2Xip^ = 6Xip J G„{k) , A/2 ^ _^ Q^^2 + 6xJ G„[k) , (22) 

fc k 

Z^ = l + \ j G„{k) , Ml=Q. (23) 

k 

Here, we have used Eq. (jA5[) from Appendix |X] and the fact that = to eliminate the pion tadpole term in Eq. 
(|17p. The numerical results for the masses, the condensate, and the pressure as functions of temperature are shown 
in the lower row of Fig. [TJ Note that Af^ = indicates that the Goldstone theorem is fulfilled at each T below the 
chiral phase transition. This is a property which does not hold in Cartesian coordinates, see for instance the lower left 
panel of Fig. [TJ Unfortunately, Af^ = also above the transition, indicating that the chiral partners do not become 
degenerate in the chirally restored phase where (p = Q. 

The effective potential (relative to its value at its global minimum) for different temperatures is shown in Fig. 
For temperatures in the transition region, one clearly observes the features of a first-order phase transition, i.e., three 
minima separated by two maxima. 



C. H ^Q,v = Q 

The chiral limit, i.e., i/ — > 0, does not exist, since for arbitrarily small, but nonzero values of iJ, the effective 
potential has no extremum near the origin Lp = 0. We have confirmed this numerically by taking successively smaller 
values of H. The reason for this non-analyticity of the chiral limit are the inverse powers of Lp in Eqs. ((T5|) " (|T8l) . 
When setting H = from the beginning, this problem does not exist because the singular terms are eliminated right 
away. Therefore, we conclude that the system cannot describe the chirally restored phase and is limited to the chirally 
broken phase. Again, polar coordinates have demonstrated a limited range of applicability. Interestingly, this problem 
does not exist in Cartesian coordinates where taking the limit iJ — > is not problematic in the chirally restored phase. 
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D. H^O,v^O 



As argued in Refs. [9|,|22l, the non-analyticity at the origin ip = encountered in the previous subsections can be 
ehminated by introducing a non- vanishing vahie for the parameter v, i.e., shifting the vacuum expectation value to 
nonzero values ip = v. Now we have to consider the full Lagrangian (|10p . In the double-bubble approximation, the 
following equations are obtained: 



= -m^((p -v) + 2X{ip -vy -H + 6X{ip ^ v) J G^(fc) 

k 



m V — H f V V 

___ + A(^-.) ^- + - 



G^k) 



k^G^ik) 



Av^X 

tp^ 



3vX 



Gaiq) 



vX 

~ V 

G.(fc) 



GAk) 



(24) 



m2 = 6A(^ - v)^ - m2 + 6A J G,{k) J k^G^{k) + 2A ^ (^3 - 2 j G^{k) , 



(25) 



Mi 



H — m V „ , V 



m V , V 



V Iv 



+ 2X-{p - vY - 2A- 2- - 3 / G„{k) 



ip \ p 



2^3 +A-j(^--)^-6A^ 



GAk) , 



(26) 



p^ 



GAk) 



(27) 



In the limit ?; — >■ oo the system of Eqs. (P^ - (P7|) reduces to the system of equations for Cartesian coordinates 
obtained in Ref. 12| [see Eqs. (28a, b), (30a) for = 2 with the replacement p ^ p — v]. This is easily explained by 
the fact that the radial as well as the angular variable are defined relative to the origin so that for large (j) we have 



0-0 



and cos 5 — 1- Hence, in the limit v oo polar coordinates become Cartesian: 



)COS — 



TT 

a COS — 



) sm — 



Solving the system of Eqs. (|24)) -p7 | in the case H ^ yields the condensate ip — v shown in the middle panel of 
the upper row of Fig. [T] (dotted and dash-dotted lines). For values v < 0.73/,r (for instance v = 0.5/^, dotted line) 
the manner in which the condensate p{T) is multi- valued reminds of the behavior for a first-order phase transition. 
However, for this to be the case, there would have to be a third solution for p — v near the origin. This solution does 
not exist and therefore there is neither a first-order transition nor a restored phase. Above v > 0.73/7r (dash-dotted 
line) the behavior of <p{T) smoothly approaches the known Cartesian behavior, as expected. 



E. The limit rricr oo 

In this section we study the properties of the nonlinear sigma model in polar coordinates at nonzero temperature. 
This investigation is also interesting in comparison to chiral perturbation theory at nonzero temperature, see also Ref. 



271 and refs. therein. 
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The nonUnear sigma model is obtained by taking the limit mo 
Eqs. (1201) and into Eq. we obtain 



oo and keeping the ratio A/m^ fixed. Inserting 



Mi 



3< 



GAk) 



1 



i^^^Sl + i G,{k) 



(28) 



For riia — > oo, the last term dominates. Therefore, if this term is positive definite, also A/^ — > oo. We have to demand 
that the last term is positive definite in order to have real- valued solutions for Since J^, G^ik) — > for — > oo, 
the positive-definiteness of the last term requires that 



(29) 



This means that the range (p < Jtv/VS cannot be described, neither for H 7^ nor for — > (i.e., — 0). Thus, 
the model is not applicable in the chirally restored phase where we expect (/? ~ 0. This result shows that it is not 
possible to describe the phase transition in polar coordinates (in the framework of the CJT formalism) when taking 
the limit to the nonlinear sigma model. 

F. Dependence on the regularization procedure 

It is important to verify how the results change when a different regularization scheme is employed. In this section 
we regularize the divergent vacuum contributions of the tadpole diagrams by applying the counter-term scheme. The 
vacuum contributions of the thermal tadpole integral 



d^k 



1 



is not neglected as in the trivial regularization scheme but is first rewritten using the residue theorem, 

d^fc 1 



(27r)^ A:2 + M2 ' 

This integral is regularized by introducing a renormalization scale ^ 



12 



d'^k 
(^ 



1 



1 



M^-M^ 

k^ + AP k^ + ^i'^ {k^ + fi^Y 



(47r)2 



Af^ln 



AP 



A'P 



(30) 



(31) 



(32) 



Figure [5] shows the temperature dependence of the condensate and the masses in the case H ^ 0, v = 0, applying 
the trivial regularization and the counter-term regularization prescription. We set fi = m„ to satisfy the constraint 
Z{T = 0) = 1 for the vacuum value of the pion wave-function renormalization factor. 

As one can see, inclusion of the counter terms does not qualitatively change the results: it is not possible to avoid 
the pion from becoming tachyonic at high temperatures. This holds also when using the counter-term regularization in 
the chiral limit. Goldstone's theorem is fulfilled since the pions are massless in the symmetry broken phase. However, 
as with the trivial regularization scheme, no chirally restored phase exists for this system when the limit iJ — > is 
taken. 

IV. THE INSTRUCTIVE EXAMPLE OF THE FREE LAGRANGIAN 



In order to explain the problems related to polar coordinates, in this section we consider a free Lagrangian in 
Cartesian coordinates: 



ree-cart = 7^^^J.$■ d'^ (j) 7^- $ ■ (j> , 



(33) 
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T [MeV] 
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FIG. 5: The condensate (left panel) and the meson masses (right panel) as functions of temperature in the case H ^ 0, v = 
in the trivial regularization (solid line) and the counter-term regularization scheme with /x = ma (dashed line). 



where > 0. The pressure can be calculated exactly: 



P free— cart 

with D^^^{k, (p) = + ml .^ = -P + . Explicitly 

V free— cart 



k 



— / dkk^ In ( 1 - . 

■jr2 



(34) 



Let us now perform the transformation to polar coordinates. The Lagrangian 

1 



(35) 



2 02-^'-- - 2 

follows directly from Eq. ([5]) by setting A = and replacing — > —m?. Note that, although the original Cartesian 
Lagrangian (j33|) is that of a noninteracting system, the transformation to polar coordinates introduces a momentum- 
dependent four-particle interaction. 

The gap equations are given by = m? , = oo, Af^ = . Because the tadpole k'^GTt{k) vanishes for AI-,^ = 0, 
see Eq. (|A5p . the pressure in polar coordinates (and in double-bubble approximation) has no contribution from the 
interaction term ^^d^T^d^TT and is given by 



Ppoiar = - I dkk^ In ( 1 - e 



TV 
90 



(36) 



The pressure p free-polar represents the sum of one particle with mass m and one particle with mass m — 0. This 
latter contribution is, however, not correct and overestimates the exact result of Eq. p4p . This simple example shows 
once more that care is needed when polar coordinates (and a certain many-body approximation) are used to study 
properties of a system at nonzero T. 



AN ALTERNATIVE WAY TO INTRODUCE POLAR COORDINATES 



The transformation 



As we have discussed in Sec. II B, the Jacobian associated with the transformation to polar coordinates introduced 
in Eq. ^ is not unity. In this section we present an alternative polar representation (-0, tt) for the 0(2) linear a-model, 
which is defined as: 



= a/^^V- cos y2^sin ^ 



(37) 
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In this case the associated Jacobian remains unity. Here the massive scalar field acquiring a nonvanishing vacuum 
expectation value in case of spontaneously broken symmetry is represented by -0. 

The Lagrangian ^ expressed in terms of the fields tfj and tt introduced in Eq. ([37| reads 



5H 



H 



19203 



+ higher-order terms. 



H 
320^ 



2^2 



(38) 



where the shift ip ^ ip + cf) has also been performed. 

The corresponding inverse tree-level propagators and tree-level masses are given by: 



™2 



-k 



2 I 2 ■ / 



H 



■i, - 16A0 + — , m; 



H 

20 



Applying the CJT formalism we derive the effective potential within the double-bubble approximation: 



(39) 



l/e//(0;G,,Gv,) -2m2 02 + 
1 



2i/0 + - / [InG-^(fc) + Z?-^(A:;0)G^(fc) - 1] 



+ ^ I [lnG-\k)+D-'{k:(b)G^{k)-l] + V,MG^,G^) 



(40) 
(41) 



^2(0;G,,G^) = -:^ [ k^G^k) fG^{l)+ 



202 

H 

6403 



G^(/c) 





^ H r 


/ G,(fc) 

.Jk 


3203 



6403 
Gv,(fc) / G,(fc) . 



(42) 



Finally, the stationary conditions for the effective potential give the following gap equations for the temperature- 
dependent masses, the condensate and the wave-function renormalization for the tjj-ficld. 



2H = im^ip + 32Xlp^ + WXtp - 



H 



GAk)-^ j G^{k) + ^ [ k'G,p{k) fG^,{l) 



A5H 

64(^4 



Gv,(fc) 



3H 

64(^4 



G,(fc) 



/ G„(fc) / G.(fc) , 



(43) 



Mj 



H 



m: = — 



H H 

2^ ^ 16(^3 

i7 



15 / G^(fc)- / G,(fc) 

Jfe Jk 



2tp 16(^3 



I G^{k) + I G^{k) , Zl = l + \ f G^{k) 

Jk Jk \ f Jk 



(44) 



(45) 



B. Results 

In this subsection we present the numerical results for this alternative polar representation. Figure [6] shows the 
condensate and the masses as function of T using the trivial and the counter-term regularization schemes in the case 
of explicit chiral symmetry breaking. In Fig. [7] the same quantities are shown in the chiral limit. In the case of 
explicitly broken symmetry there is a crossover transition and one observes no degeneration of the chiral partners at 
high T. Just as in Sec. Ill, for T > 300 MeV the validity of the model breaks down and the numerical results are no 
longer reliable. Moreover, chiral restoration is approached very slowly, also when the chiral limit is taken. 
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FIG. 6: The temperature-dependent condensate (left panel) and the meson masses (right panel) for the alternative polar 
representation in the case of explicitly broken symmetry, in the trivial regularization (solid lines) and in the counter-term 
regularization scheme (dashed lines) with fi = m^. 




FIG. 7: The temperature-dependent condensate and ^-mass for the alternative polar representation in the chiral limit, in the 
trivial regularization (solid line) and in the counter-term regularization scheme (dashed line) with fi — m^. 



Summarizing the results the alternative polar representation does not offer a better description of the chiral phase 
transition (at least not within the chosen many-body approximation). The reason for this is that, vifhen constructing 
this alternative polar representation, one has to evaluate several terms of the form /(I -I- x) with x = ip/4> and 
X = 7r/(/). Increasing the temperature the condensate starts melting while the fluctuation of the fields become larger 
and the calculation becomes no longer reliable for higher temperatures. Especially when approaching the critical 
temperature all the terms of the form x/(j), with x ~ TT,ip, become problematic, limiting the validity of the model to 
lower temperatures. 



VI. CONCLUSIONS 



In this paper, we have studied the 0{2) model in polar coordinates at nonzero temperature. After having clarified 
some issues related to the transformation from Cartesian to polar coordinates in the functional integral representation 
of the partition function, we have computed the latter in the C JT formalism in double-bubble approximation. We have 
studied in detail the cases where the chiral symmetry is explicitly broken and where the explicit symmetry breaking 
parameter is set to zero. We have distinguished the latter case from the chiral limit, where the explicit symmetry 
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breaking parameter is smoothly sent to zero. We have found that this limit does not exist in the strict mathematical 
sense, due to non-analytic terms 

\-7^dn,Trd^iT and Ha cos — 
2 (j)^ (p 

appearing in the Lagrangian for polar coordinates. For all cases except when the explicit symmetry breaking parameter 
is exactly zero, we have found that, when approaching — > 0, the ensuing divergences are sufficiently severe to 
invalidate the approach above a certain maximal temperature T^ax- Above T^ax no physical solutions exist. The 
same results hold in the nonlinear limit mo- — > c». 

We have also investigated the possibility of shifting the potential by an amount v along the (/)i-direction in order to 
circumvent the divergences resulting from the non-analytic terms in the Lagrangian. Variation of v allows to change 
smoothly from polar coordinates to Cartesian coordinates (corresponding to u — > oo). Our conclusions about the 
suitability of polar coordinates remain unchanged, however, T^ax increases with u, so that the range of applicability 
of polar coordinates is extended. Morevocr, the larger v, the better the agreement with the Cartesian results. 

We have also introduced and studied an alternative representation for polar coordinates (with unit Jacobian in the 
functional integral representation of the partition function). Also in this case, the above mentioned problems persist. 

Our most important result is that, both in the chiral limit and in the case of explicitly broken symmetry, the chiral 
partners do not become degenerate in mass at high T, not even when approaching T,nax where the order parameter 
has already decreased by a substantial amont. In general, the sigma particle becomes more massive while the pion 
mass decreases or remains zero (in the absence of explicit symmetry breaking). Above T^ax the pion even becomes 
tachyonic. The absence of degeneracy of the chiral partners means that an important indication for the restoration 
of chiral symmetry is missing when using polar coordinates. We conclude that the use of angular variables is not well 
suited for the study of the chiral phase transition. 

Note that we used either trivial regularization or the counter-term renormalization scheme when treating the 
divergent vacuum contributions to the thermal integrals. Our results are qualitatively independent from the chosen 
regularization scheme. Due to consistency with the extension of the range of functional integration in polar coordinates, 
one should perform the renormalization procedure in dimensional regularization, which remains a task for the future. 

A possible extension of the present study would be to use four-dimensional polar coordinates, corresponding to the 
0(4) model. The polar 0(4) model has the advantage that three degrees of freedom can be identified with the three 
pions, TT*^ and tt^, and the remaining one with their chiral partner, the sigma particle. The number of angular degrees 
of freedom could affect the behavior of the equations in the limit iy9 — )■ 0. However, we believe that this generalization 
will not fix the problem encountered in the polar version of the 0(2) model, since the divergences pointed out above 
are general. 

Finally, although on a conceptual level each representation is equivalent, on a practical level the use of Cartesian 
coordinates is favorable to study thermodynamical properties of systems described by the 0{N) model. In connection 
to QCD, one should extend the model by incorporating all the relevant low-energy mesons: besides scalar and 
pseudoscalar particles, also vector and axial- vector degrees of freedom should be included in an enlarged U{N) x U (N) 
symmetry for a more realistic treatment of properties of QCD at nonzero T 
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Appendix A: Thermal integrals 



We list the standard thermal integrals which were used for numerical calculations. Carrying out the Matsubara 
summation of the thermal tadpole integral gives 



1 



cxp 



(Al) 
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consisting of a finite contribution, 



dk 



and a divergent vacuum contribution, 



m = 



d^k 



d^k 1 



For the explicit calculation we need the following integrals: 

j G, = Q^{M„) + Qt{M„ 



Z 



Z 



k'G^ 



Z4 



Z 



Z 



For the effective potential we need in addition: 



i-k' + ml)- 



1 



-fc2 + Af2 -fc2 + M2 



= ^{ml - AC) I ^^^^ = ^iml - Ml) [Q^{AU) + Qt{M,)] 

k 



{-k^ + ml) 



\{<-Ml)^, 



Ah 
Z 



k 

Qt 



Z 



1 -Z2fc2+M2- 

-Z'^k^ + Afi -Z^k-2+AI^_ 

Qt 



2^ ' Z^ 



AI^ 
Z 



Z 



\ I \nG-' = ^R^{AU) + ^Rt{AI,) 



-Rt 



AI^ 
Z 



where 





which in the case M = and Z simplifies to 



i?T(0) 



'4^2 



^5 



(A2) 

(A3) 
(A4) 
(A5) 



(A6) 



(A7) 

(A8) 
(A9) 

(AlO) 
(All) 



Except for the discussion in Sees . IIII Fl and IVl we neglected the contributions from renormalization, i.e., = as 
well as R^ = 0. This approximation scheme is called the trivial regularization. 
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